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Abstract 



Integral relations and transformation rules are used to obtain, out of an asymptotic 
solution, a new group of four pairs of solutions to the double-confluent Heun equa- 
tion. Each pair presents the same series coefficients but has solutions convergent 
("^ ] in different regions of the complex plane. Integral relations are also established 

CN ' between solutions given by series of Coulomb wave functions. The Whittaker-Hill 

equation and another equation are studied as particular cases of both the double- 
, confluent and the single-confluent Heun equations. Finally, applications for the 

I Schrodinger equation with certain potentials are discussed, mainly for quasi-exactly 

Oh' solvable potentials which lead to the above special equations. 

> 



1. Preliminary remarks 



Here we deal with two groups of series solutions to the double- confluent Heun equation 
X ■ (DCHE). The first group, constituted by four pairs of solutions, is generated from an 

■ asymptotic expansion by means of integral relations and transformation rules, and the 

second group is given by pairs of solutions in series of Coulomb wave functions, already 
derived in [1]. For the latter, we show that an integral relation also exists between the 
members of each pair, and we provide additional properties for the solutions. After this, 
we analyse two differential equations which are special cases of both the DCHE and the 
single-confluent Heun equation and, finally, we use some results to solve the Schrodinger 
equation for certain potentials. Before proceeding, we set down some conventions, present 
the procedures used to obtain the solutions and outline the structure of the paper. 
For the DCHE we adopt the form 

z^^ + {B, + B2z)^ + [B,-2r]u;z + ujh')U = 0, ^ 0, oo ^ 0) (1) 

where z = and z = oo are irregular singularities. Bi = and/or uj = are excluded 
because, in these cases, the equation reduces to a confiuent hypergeometric equation or 
to an equation with constant coefficients. Thus, if 5i = and u ^ 0, the substitutions 



y 



-2iujz, U{z) = e-y/^^fiy), - (1 - Es)^ + ^3 = (2a) 
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give the confluent hypergeometric equation 



,_^ + [(2„ + B,)-rf-^-(,, + „ + ^)/ = (2b) 

and, if 5i 7^ and a; = 0, 

y^Bjz, U{z)^y^g{y), /^^ - (S2 - + S3 = ^ (3a) 
^0 + [(2/5 + 2-S,)-y]|-/?, = O. (3b) 

If Si = a; = 0, we find an equation with constant coefficients by taking z — expy. 

Equation (1) can be obtained from the generahzed spheroidal wave equation (GSWE) 
[2], which is also known as single-confluent Heun equation or, simply, confluent Heun 
equation [3]. Actually the GSWE in Leaver's form reads [4] 

(PU dU r n 

z{z + {Bi + B2Z)— + [S3 - 2r]uj{z - Zo) + uj^z{z - Zo)\ = 0, {uj ^ (4) 

where zq^ Bi, t] and u are constants, z = and z = zq arc regular singular points while 
z = cxD is an irregular singularity. For Zq = the preceding equation gives the DCHE and, 
when a limiting solution also exits, this one can be taken as a starting-point to generate a 
group of solutions to the DCHE. This procedure was applied by Leaver to find expansions 
in series of Coulomb wave functions for the DCHE, although these can as well be found 
otherwise [5]. Note that, although the DCHE was called 'confluent GSWE' in [1, 4], it is 
not derived from the GSWE by a process of confluence (see [6], chapter 4). 

The transformation rules aforementioned result from variable substitutions that con- 
vert the DCHE into another version of itself. Thus, for a given solution S{z), 

S{z) :=U{B,,B2,Bs;u;,rj-z), (5) 

(':=' means 'equal by definition') we shall generate new solutions by using the following 
rules [1,3,6-8] — ri, r2, rs — 

nS{z) = e''^"+^i/(2")z-'''-^^/'C/(S;,S;,S;;a;',r/';i?), (6) 

r2S{z) = e^i/^^2-^^C/(-Si, 4 - S2, S3 + 2 - B^; u, 77; z), (7) 

r3Siz)^UiBi,B2,B3;-uj,-r);z), (8) 
where, on the right-hand side of the first relation, we have 

B[^u;B,, S; = 2 + 2ir7, s; = S3 - + zr/) - tr] - 1 

u; = 1, z?7 = — - 1, 1? = — . (9) 

In the third rule it is assumed that we must change the sign of (?7,a;) only where these 
quantities appear explicitly in S{z), keeping the expressions for the other parameters 
unchanged, even if they depend on 77 and cu. Moreover, note that the transformation ri 
generally gives a solution having a region of convergence different from that of the initial 
solution, since it involves an inversion of the independent variable. For brevity, we use 
only ri and r2, although must be employed to furnish the full group of solutions. 
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An integral relation can also be used to generate a new solution from a previously 
known one. It provides a pair of solutions which has essentially the same series coefficients 
(that is, coefficients which differ at most by a constant not depending on n) but put some 
restrictions on the parameters and arguments of the solutions. The transformation rules 
applied to that pair afford new ones, having again the same coefficients and solutions 
connected to each other by integral relations. 

We denote by the series solutions which converge for \z\ > 0, and by f/° those which 
converge for \z\ < oo, the subscript i indicating the pair to which the solution belongs. 
If there is a phase parameter u, it appears as a subscript as well, let us say, and 
For a specified pair, C/° will result from via an integral transformation. Solutions 
with a phase parameter are given by two-sided series in which the summation index n 
runs from — oo to oo, whereas solutions without that parameter are given by one-sided 
series (n > 0). Under certain conditions, the latter become finite-series solutions which 
are called quasi-polynomial solutions, Heun polynomials or quasi-algebraic solutions. 

In the next section we find the kernels to the integral relations used in sections 3 
and 4. In section 3, we use integral relations and transformations rules to derive four 
pairs of sohitions from an asymptotic expansion in the vicinity of ^ = oo. In each pair, 
one solution is given as an ascending or descending power series of z and the other as a 
series of irregular confiuent hypcrgcomctric functions which, for terminating series, may 
be written in terms of generalized Laguerre polynomials. 

In section 4 we determine the integral relations between expansions in Coulomb wave 
functions. We use only expansions in series of irregular hypergeometric functions, and 
find that, in each pair without a phase parameter, one solution may again be expressed 
as a generalized Laguerre polynomial. In addition, wc verify that finite-series solutions 
occur under the same conditions valid for the corresponding solutions in section 3. 

In section 5 we obtain normal forms for the DCHE — in which there is no first 
derivative term — and examine the two differential equations which are particular cases 
of both the DCHE and the GSWE, namely 



dv? 
d'^W2 
dv? 



+ [Oo + 9i cosh(fi;M) + 02 cosh(2fi;M)] Wi = 0, (10a) 
$0 + 9i smh{Ku) + 02 cosh{2Ku)] W2 = 0, (10b) 



where the first equation represents the Whittaker-Hill equation (WHE) or the modified 
WHE depending on whether nu is pure imaginary or real, respectively. We also discuss 
the solutions to the Schrodinger equation for some potentials which give rise to DCHEs, 
and specially to these particular cases. Section 6 is devoted to final comments while, in the 
Appendix, integrals used to establish integral relations are written in terms of irregular 
confluent hypergeometric functions rather than in terms of Whittaker functions. 



2. Kernels for integral relations 

Several kernels are possible for integral transformations of the DCHE [5, 6], but we only 
regard those that will be useful in the subsequent sections. We follow a procedure similar 
to the one employed by Schmidt and Wolf [5], adapted to the form we have chosen for the 
DCHE. Thus, if U{z) is a known solution of equation (1), we seek a new solution U{z) 
given by the integral relation 

U{z) = r K{z,t)U{t)dt, (11) 
Jti 
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where the kernel K{z, t) is determined from [9] 

L,{K{z,t)} = Lt{K{z,t)}, 
being the operator and its adjoint given by 

:= z% + [B, + B2Z] £ + {uo'^z^ - 2ujr]z) , 

Lt ■^t% + [-B, + (4 - B2)t] I + {coH'' - 2u;rjt + 2-B2). 
In terms of L^, equation (1) reads 

[L, + Bs]Uiz)^0, 



(12a) 



(12b) 



(13) 



where is now understood as an ordinary differential operator. In equation (11) we have 
chosen the contour of integration as the line joining ti and t2, but we assume that these 
endpoints depend on z and, consequently, we have to use the formula 



d 

dz Jti 



,-,/ N , [^■^ dFiz.t) , ^, ,dt2 ^, .dti 



(14) 



in order to derive the conditions under which U{z) is solution of the DCHE. Hence, 
applying to integral (11) and using equations (12a) and (14), we find 



LAUiz)} = r Lt{K{z, t)}Uit)dt + Qiz, h) - Qiz, h), 
Jti 



(15a) 



where {i — 1,2) 
Q{z,ti) 



U{ti)K{z,ti) + 

dK{z,U) fduV ^ ^ dK{z,U) dti 
dti \ dz I dz dz 



/duy duiu) 

\ dz I dti 



K{z,ti). (15b) 



The notation Lt{K{z,t)}U{t) in equation (15a) means that the operator Lt acts uniquely 
on the object inside braces. Next, we integrate equation (15a) by parts or, equivalently, 
by using the identity 

Lt{K{z,t)}U{t) - K{z,t)Lt{U{t)} = ^P{z,t), 
where P{z,t) is given by 

P{z,t) -.= 1^ 

Then we find 



dt 



dt 



+ [{2-B2)t-Bi]U{t)K{z,t). (16) 



L,{Uiz)} 



t2 



K{z,t)Lt{U{t)} + 



dP{z,t) 



dt + Q{z,t2) -Q{z,ti), 



or, using equations (13) and (11), 

(L, + B3)U{z) = P{z, t2) + Qiz, t2) - P{z, ti) - Q{z, h). 



(17) 
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Therefore, when U {z) is a solution of the DCHE, U [z] will also be a solution if the integral 
(11) exist and the right hand side of equation (17) vanishes. In fact, we require that the 
'integrated terms' vanish when t ^ U, that is. 



Q{z,U) + P{z,U)^0, 
where, due to equations (15b) and (16), 
P(z,ti) + Q{z,ti) = 

.2 



(18a) 



dti 



dz dz 



K{z,ti)U{ti) + 

dK{z,t 



(I) - 



X 



U{U). (18b) 



Below, we find two kernels for integral relations and show that, for these, the right-hand 
side of (18b) assumes a very simple form. From these kernels, two others — and the 
corresponding transformed solutions U{z) — may be obtained by the change (77,0;) 
(—77, —Lo). These kernels are necessary to set up integral relations between the solutions 
derived from the ones considered in this paper by means of rule r^. 

First kernel: Performing the substitution 

K{z, t) = e^'^(^+*)+(-BiA) z^-^'Hi{z, t) (19) 

in equation (12a) and supposing that Hi depends on z and t by the product zt, we find 

2) Hi^O, ^ := -2iujzt/Bi =^ 



(e-l)^-(f -zr^ 

Hi{z,t) = (^ - l)(B2/2)-i^-2 ^ (^-^zt - 1) 

Thus, we have a first kernel, denoted by Ki{z,t), namely. 



^-i„-2 



(20) 



Ki{z,t) = e 

which yields 

Ui{z) = e 



iuj{z+t)+{Bi/z) ^2-^2 



2iu! 



Bi 



■zt - 1 



(21) 



*2 



iijt 



2iuj 



Bi 



-zt - 1 



B2 ■ o 



We rewrite this integral in terms of ^ and integrate from ^1 to ^2, assuming that these 
new endpoints are constants to be specified later. We get 



On the other hand, by inserting 



(e-1) 



(B2/2)-ir?-2 



u 



2iujz 



(22) 



ti{z) = -Biii/{2iuz) 



(23a) 
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and K{z,t) — Ki{z,t) into the right-hand side of the expression (18b), we find that the 
first term becomes 



{B,)% B,UB2 - 2) ^ 
H : -Di 



lUJZ 



Ki{z,ti)U{ti), 



the second term vanishes because z'^ {dti/dz)^ — {ti)^ — 0, and the last term reduces to 



2iu!z'^ 

Therefore we can write 

Pi{z,ti) + Qi{z,t 



ILOZ 



Ki{z,ti)U{ti). 



+ Bi 



{Ci-l)K,{z,U)U{U 



2iujz'^ 



+ Bi 



2iijjz^ 



ILOZ H 



z 2z 



. (23b) 



Now we choose = 1 ^^'^ ^2 = 00 and then we can use one of the integrals given in the 
Appendix to integrate equation (22) for the solutions given in sections 3 and 4. 



Second kernel: Accomplishing the substitution 

K{z,t) = e"^^'+'^-^^'/'h'''-^H2{z,t) 



(24) 



in equation (12a)and supposing that H2{z,t) depends on z and t by the product zt, we 
find 



(C-l)^ + (i?7 + f )/^2 = 0, C-^'2iu;zt/B 
H2{z,t) = (C - l)-(^2/2)-i,, ^ (^^^^ _ 
Hence, the second kernel, K2{z,t), is 



B2 ■ 



(25) 



K2{z,t) = e-(-+*)-(^i/*) t^--' (^zt - 1 

V Bi 



-IT) 



(26) 



which implies 



U2{z) = e"^' dt 



or, in terms of C, 

UJz) = e^'^^Z^-^^ dC L«iC/(2^)-2WC^B,-2 _ ^y{B,/2)-ir, ^ f B^C 

JCi [ \2iujz^ 

This time we have 



(27) 



ti = B^Q/{2iuz) 



(28a) 
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in the condition (18a), and equation (18b) can be rewritten as 



P2{z,ti)+Q2{z,ti) 



2iujz 



B2-2 



{Ci-l)K2{z,U)U{U)-- 

Bid 2,100 z 



{Bi)% 
2iuiz'^ 



-Bi 



X 



lUJZ + 



2z 



(28b) 



We choose Ci = 1 ^^nd (^2 = 00 and, then, we can once more use one of the integrals of the 
Appendix to integrate equation (27). 

The approach we have used in this section is similar to the one employed by Schmidt 
and Wolf [5], in which we have regarded ^ or C as integration variables instead of t. This, 
in turn, implies that ti and t2 in the the integral (11) are functions of z, as far as the 
integration endpoints in the variables ^ or C are taken as constants throughout sections 
3-4. As said before, we shall find U{t) = U°^{t) and lJ{z) = U^{z) in the integral relations. 
Then, using the first identification together with the validity conditions for the integrals, 
it will be easy to check that the right-hand sides of equations (23b) and (28b) vanish for 
the solutions given in sections 3 and 4. 



3. Solutions derived from an asymptotic expansion 

The simpler quasi-polynomial solutions to the DCHE are obtained from the asymptotic 
expansions in the vicinity of the singular points or cxo [5] . In this section, we begin with 
an asymptotic representation for 2; —> 00 and use the integral relations and transformation 
formulae to form four pairs of solutions. In each pair, one solution is given in terms of 
a series of ascending or descending powers of z, and the other in terms of a series of 
irregular confluent hypergeometric functions. We flnd the conditions for obtaining Heun 
polynomials and verify that, in this case, the hypergeometric functions degenerate to 
generalized Laguerre polynomials. 

The starting-point solution is given by 



(29a) 



n=0 



where the recurrence relations for the coefficients b!^^ are, in abbreviated notation. 



aoh + Pobo = 0, anbn+l + Pnbn + Tn^n-l = (n > 1), 



(29b) 



in which 



B2 

2 



+ trij(^l + irj-^y (29c 



Pn^n{n + 1 + 2iri) + iuBi + B3 + 
7n = 2iujBi [n + iri+ {B2/2) - 1] . 

These relations yield a characteristic equation in terms of the infinite continued fraction 

ao7i Q;i72 "273 



^° f32- Ps- ■ 

To obtain the foregoing recurrence relations, we perform the substitutions 

U^{z) = e'^'z-'''-^^'/^'^Y{y), y = {-2icuz)-\ 



(30) 
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and find 



2^ 



+ 



1 + (2 + 2ir))y + 2iujBiy^ 

B. 



dY 
dy 



+ 

So 



C + 2iuBi irj + 



Bn 



y = o, 



C := ^a;Si + ^3 + + zr;) + ir? - ^) . 



Then, inserting Y{y) = J2'^=ob^n ^y^ into this equation and proceeding in the usual form, 
we obtain the previous relations. 

According to the theory of ordinary differential equations in the complex domain (see 
[10], chapter 7), the solution U^{z) is unique (one- valued) within the sector 



— ^ < &rg{-2iujz) < —. 



(31) 



However we still have to show that it converges for |2;| > 0. To accomplish this, we divide 
the recurrence relations (29b) by and take the limit when n — > 00. This gives 



hm = or hm 



n 



n— >oo 



—n. 



(32) 



These limits may as well be derived by using a Perron-Kreuser theorem for difference 
equations [11]. In order to satisfy characteristic equation (30), we have to choose the first 
limit (minimal solution). Thence, 



\Sl^{-2iu:z) 



lim 



-n-l 



Si 



}&^{-2iu:z) 



-n nz 



and, therefore, U^{z) converges for any \z\ > 0. 

On the other hand, to get a solution U^{z) convergent in the neighborhood of 2; = 
we insert U^{t) into the integral relation (22). This gives 



n=0 



being 



\ R /. /Si\i+"+2*'' / B2 Bi 

-iV-n e'^'^' i-j] [/ f n + i77 + ^, n + 2 + 2ir), -j- 



B2 
2 



where we have used integral (Al) which is valid if 

3? [(S2/2) - ir/ - 1] > 0, ^{Bi/z)>0. 
Thus, the sought solution is given, apart from a multiplicative factor, by 

00 ^ rj R \ 

n=0 ^ 2 ^ ^ 



(33) 



(34) 



where U{a,b,y) denotes the irregular confluent hypergeometric function [12]. Note more- 
over that, by inserting into equation (23b), we have 
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Pi{z,ti) + Qi{z,ti) oc 



X exp 



iuz - — (Ci - 1) 

z 



2iLiJz'^ 

n=0 \ ^ 



,2+i7j-(S2/2)^7i'?-(S2/2)^^^ _ ^^(B2/2)-i7]-l 



The series on the right-hand side converges at .^i = 1 and at ^2 = 00. Then, the condition 
(33) assures that, for ^1 = 1, the second member goes to zero since (^1 — l)(^2/2)-i?7-i _^ q. 
for ^2 = 00, the second member also vanishes because exp [— i?i(^2 ~ l)/-^] —>^0. Argu- 
ments similar to these may be repeated for the other pairs of solutions written below. 
To obtain the behaviour of U^{z) when 2; — > 0, we use the relation [13] 

U{a,b,y)-^y-''[l + 0{\y\-\ -37r/2 < arg|/ < 37r/2, (b| ^ 00), (35) 
thereof we find 

-|- ) , within the sector - — < arg ( — ) < — . (36) 
BiJ 2 \ z J 2 

However, to show that the series in Ui{z) converges for \z\ < 00, we must consider the 
behaviour of U{a,b,z) when b — > 00, while b — a and z remain bounded [13]. In this 
manner we get 

lim ■^"''"^ — ^ ^ rr , „•„ , -^2 _ , o , o„„ 



= fn:=U (n + i7] + ^,n + 2 + 2tr],^y 



fn 

Then, combining this expression with the first limit given in (32), we have 

lim ^Si^ = (37) 

Therefore, Ui{z) converges in any finite region of the complex plane. In case of finite series, 
the ratio test becomes meaningless and the convergence must be decided by inspection. 

Starting from the first pair of solutions, we generate three others by the transforma- 
tion rules ri and r2 and, in each pair, the solutions are also connected by an integral 
transformation. Below, we collect up the four pairs of solutions (^7°°, U^) {i = 1, 2, 3, 4), 
the validity conditions for the integral transformations, the asymptotic behaviour of each 
solution, and the sufficient condition to obtain quasi-polynomial solutions. For the latter 
solutions, the the functions U (a, b, z) degenerate to an generalized Laguerre polynomials 
because the parameter a becomes a negative integer — /, and thus we have [12] 

[/(-/, 1 + a, y) = (-l)'Z!Lr(y). (38) 

The condition for quasi-polynomial solutions results from the fact that a series with three- 
term recurrence relations such as (29b) becomes a finite series with 0<n<A'" — lif 
7„ = for some n — N ^positive integer [14]. 

First pair: {U^,U^). 

U^(z) = £^'^^^-^''-(■^2/2) ^^^^ b^^\-2iu;zy, 

(39a) 

U^(z) = e^^^^-'^-(^2/2) ^oc^^ i){i)u (^n + iri + ^,n + 2 + 2ir], . 
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= n (n + 1 + 2ir]) + iuBi + B^+[^ + tr]j i^l + tr]-^j, (39b) 
7^^) = 2iujBi (n + iri + ^-iy 

Integral relation (22) : 3?(^^ - «?7 - l) > and > 0. (39c) 

Finite series: (S2/2) +i77 = l- A^ ^ 0<n<A^-l. (39d) 

To find the integral relation we have used equation (Al). The asymptotic behaviours 
are 

lim U^{z) ~ e''^'z-'''-^^^/^\ -— < arg {-2iLUz) < — ; (40a) 

lim [/{>(.) ~ 1, -y < "'Kf ) < f • 

Second pazr : oc rsC/f" and C/2° « '^2^/1°. 
[/2°°(^) = ei'^^+(-BiA)z-i'?-{-B2/2) ofeW(_2ia;^)-'^, 

[/O(^) = ei'^^+(SiA)^-i'?-(S2/2) ^^^^ i){^)u (n + 2 + i7] - f,n + 2 + 2ir], 



(41a) 



pj^^ = n (n + 1 + 2i7^) - iuBi + S3 + ^ + + ~ t) ' ^^^^^ 
7i') - -2^a;i?i (n + 1 + ^r^ - ^)) . 

Integral relation (27): 3?(^^ + i?? - l) < and ^(^^ < 0- (41c) 

Finite series, if {B2/2) - irj = 1 + N =^ 0<n<N-l. (41d) 
The integral relation is found by using equation (Al). The asymptotic expressions are 

Jim U^(z) ~ e^'^^z-^^-^^^/'), -y < arg (-2iu;z) < (42a) 

lim t/2°(^) ^ < arg (-^) < ^. (42b) 



Third pair : oc rit/° and oc riU^ . 

U^iz) = e-^ E~ b^n^U (n + i77 + f , n + S2, -2iu;z) , 



(43a) 



ai^) = n + 1, 

/5^3) = n (n + ^2 - 1) + icu^i + S3, (43b) 
7^^) = 2ia;Si (n + ir] + ^ -1^ . 

Integral relation (22): 3?(^^ - iry - 1^ > and > 0. (43c) 

Finite series, if (S2/2) + iri = I - N ^ 0<n<N-l. (43d) 
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To find the integral relation we use the equation (A2). The asymptotic behaviours are 
the same as in the first pair, that is, 

jim U^{z) - e*--^-*'?-(S2/2)^ _^ < arg {-2iujz) < y ; (44a) 

lim[/3°(.)~l, -T<^^s(f)<f- ^''^^ 

Fourth pair : oc r2U^ and oc r2U2- 

U^{z) = e^'^^+(-^i/^)^2-B2 ^^^^ h(A)u(n + 2 + ir)-f,n + A-B2, -2iuz) , 



(45a) 



= n + 1, 

^ ^ + 3 _ 5^) + 2 - - ^2 + B^, (45b) 
7^^^) = -2ia;Si (^n + 1 + ir? - ^) . 

Integral relation (27): 3fj(^^ + ir]-l^ < and < 0. (45c) 

Finite series, if (S2/2) -ir]^l + N =^ 0<n<N-l. (45d) 

To find the integral transformation, we must use relation (A2) again. The asymptotic 
behaviours are the same as in the second pair, namely, 

jim U^{z) ~ e^'z-''^-^^'/^\ -y < arg {-2iujz) < y , (46a) 
lim t/°(.) ^ e^^/^.^-^^ < arg < ^. (46b) 

For 2; — > 0, we have found the two asymptotic behaviours we could expect from the 
theory of differential equations. However, for z 00, we have only one of the expected 
expressions. This occurs because we have regarded only one half of the solutions. In 
effect, if we apply rule rs to the preceding solutions, we get four new solutions U°°{z) 
(i = 5, • • • , 8) for which 

jim U^{z) ~ e-^'^^^^''-^^^/"), < arg {2iujz) < y . 

Note as well that U^{z) and U^iz) can be reexpressed in terms of confluent hyperge- 
ometric functions. Indeed, if we use integral (A2) to derive these solutions from U^{z) 
and U^{z), we find 

Uliz) oc e-^z-^^-(^^/2) f^b^^)ufn + irj + ^,n + l + irj + ^, ^) , 

n=0 ^ 2 2 Z J 

Ul{z) oc e--+(^i/-)^-^-(«^/2) g 6(4)c/ + 2 + i77 - n + 3 + i77 - ^, , 

n=0 ^ 2 2 Z J 

which are consistent with the previous expressions due to the transformation [12] 

[/(a,6,y) = ^^-^t/(l + a-6,2-6,|/) (47) 
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followed by (38) with / = 0. Prom the above expression for U^{z), we find 

b);'U[n + irj+ ' • 

fi=0 

U^{z) = r2U^{z) oc e^"^+(^i/^)^2-S2 ^ 5(2) 



U^{z) = riU^iz) oc e'"^' (n + ir) + + 1 + ir) + -2itjz) 



00 
n=0 

/ B2 B2 \ 

U \ n + 2 + ir] n + 3 + ir] —2iujzJ , 



which are alternative forms for U^{z) and U^{z). It is useful to write the two solutions 
in each pair as series of hypergeometric functions because we can use equation (35) to 
deduce the asymptotic behaviours of both solutions as well as the respective sectors inside 
which they are one- valued. 

4. Solutions in series of Coulomb wave functions 

In this section wc establish integral relations for expansions in scries of Coulomb wave 
functions (already found in section 4 of [1]) and provide some additional properties for 
these solutions. We consider only expansions in series of irregular confluent hypergeo- 
metric functions and, thus, we use the notation {U°°, U^) where we have used (t/j, Ui) in 
[1]. In effect, the solutions in terms of U{a,h,z) afford the expected behaviour for the 
solutions when z — > 00 and — > 00. Moreover, in section 4.2 we note that, in each pair 
of solutions without a phase parameter, one solution may be expressed in series of gen- 
eralized Laguerre polynomials, and we also find that the conditions for quasi-polynomial 
solutions are the same as in the corresponding solutions of section 3. Therefore, to discard 
the expansions in regular confluent hypergeometric functions does not imply that we are 
setting aside flnite-series solutions. 

4.I. Solutions with a phase parameter 

The flrst pair below is equivalent to the solutions found by Leaver [4]. U^i,{z) can be 
derived from U^{z) by the rule ri and also by an integral transformation. The second 
pair results from the first one by means of the rule its solutions are connected to one 
another by an integral relation but not by the rule ri. 



(48a) 



First pair : 

U^^{z) = e*'^^^'^+i-(-B2/2) Y.n=^oc hn{-2iujzYU{n + u + 1 + iri,2n + 2u + 2, -2iujz), 
Ul{z) = e-^^-'^-^^^/^) E~ _oo bn (f U(n + u+f,2n + 2u + 2,f), 
with the following recurrence relations for the coefficients 6„ 

ttn&n+l + Pnbn + InK-l = 0, (48b) 



where 



_ iwB-i [n-\-v+2-{B2/2)] [n-\-v+l-ir)] 



/3„ ^ ^3 + (n + . + 1 - f ) (n + . + f ) + ^ (48c) 

_ icjgi [n+v+{B2/2)-l][n+v+irf\ 
2[n+v][n+v-{l/2)] 

Integral relation (22): 3?[(52/2) - ir; - 1) > 0, ^{Bi/z) > 0. (48d) 
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The phase parameter v may be determined from a characteristic equation given as a sum 
of two infinite continued fractions, namely, 

_ q;_i7o q;_27-i Qi-37-2 OoTi 00^0273 

" /3-2- /3-3- ■ ■ /32- /33- ■ ■ ■ ■ ^ ^ 

To show that these solutions are connected to each other by the integral relation mentioned 
above, we introduce U{t) — U^{t) into the integral (22). Then we find 



CXD 



since := —2iujzt/Bi) 

/« = e-^i^/^ - ^„+.+i_(s,/2)^ (^^ + ^ + 1 + 2n + 21/ + 2, 

= r - i77 - 1) e-^^/^[/ (^n + 1/ + ^, 2n + 2i/ + 2, , 

under the conditions written in (48d). Furthermore, by inserting U^{z) into condition 
(23b), we find {U = -Biii/{2iujz)) 



Pi{z,ti) + Qi{z,ti) oc 



cm 

2i6i;2;^ 



z 



l_^_(B2/2)^l+!/-(S2/2)^^^ _ -|^^(S2/2)-i7?-l 



X exp 



iujz -—{ii- 1) 
z 



E^nM t/(n + z/+l + i77,2n + 2j/ + 2,— 

n=0 \ ^ / ^ ^ 



As in section 3, the right-hand side of this expression vanishes for ^1 = 1 and ^2 = 00 
and, therefore, the condition (18a) is satisfied. 

Second pair : 

U^iz) = fiz)z''+^-^^^/^'^ E^=-oo b'^i-2iu;zruin + ly + 1 + ir),2n + 2iy + 2, -2iu;z), 
Ul{z) = /(^)^-^-(^^/2) E^=-oo t (-f (n + + 2 - f , 2n + 2i. + 2, -f ) , (50a) 

where the recurrence relations for 6^ are 

«n^n+l + ^n^n + ln^n-\ = 0, (50b) 



with 



' _ lUiBi [n+v+{B2/2)][n+v+l-irj] 
"n — 2[n+v+l][n+u+Ci/2)] ' 



= -ij3 - (» + + 1 - f ) (n + . + f ) - TntK I'lS' ■ W 

' _ iujBx [n+v+\-{B2 /2)] [n+t/+iT?] 
'''n ~ 2[n+i/][ra+!^-(l/2)] 

Integral relation (27): 5R[(S2/2) + ir^ - 1) < 0, 5R(Si/z) < 0. (50d) 

The characteristic equation is analogous to (49). Again, it is simple to find the integral 
relation stated above. We insert U^{t) into the right-hand side of relation (27) and get 

oo / R \ 

U,{z) oc e^-.— (^V^) Y: b: (-^) /f (.) oc Uliz), 
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since {( :— 2iujzt/Bi) 

7(2) ^ /- e^i</^ (C - (-n+u-i+{B2/2)u + ^ + 1 + 2n + 2z/ + 2, -f C) 

= r (l - i?7 - ^) e^i/^C/ (n + 1/ + 2 - ^, 2n + 21/ + 2, . 

On the other hand, by inserting into the equation (28b) we find that the condition 
(18a) is satisfied. Moreover, wc note that the vahdity conditions for the integrals do not 
involve the phase parameter and, therefore, these integral relations are also valid for the 
truncated solutions. 



Ji.-2. Solutions without phase parameter 

These come from the truncation of the solutions with phase parameter (n > 0) but, 
contrary to the solutions given in section 3, now there are three possible forms for the 
recurrence relations and for the corresponding characteristic equations (see appendix of 
[1]). For completeness, we write out these relations. In the first one we have q;_i = and, 
in the other cases, q;_i may be different from zero. 



aohi + /^o&o = 0, 

anhn+l + (inbn + InK-l = (n > 1), 

Oiobi + /3o&o = 0, 

0:162 + Pibi + [q;_i + 7i] 60 = 0, 

ttn&n+l + f^nK + InK-l = (n > 2), _ 

(y-ohi + [/3o + a-i] &o = 0, 

anK+l + PnK + lnK-1 = (n > 1), 



/3o 



/9o 



Qo7i Q:i72 0:273 
^ ^ ^ 



ao [ot-i + 71] q;i72 q;273 



/3o + o- 



ao7i 0:172 0273 
^ ^ 



(51) 



(52) 



(53) 



In each one of the the four pairs of truncated solutions, one solution can be expressed 
in terms of generalized Laguerre polynomials by the Kummer transformation (47) followed 
by equation (38). Such solutions are: U^{z), U2{z), U^{z) and U^{z). In the case of 
finite-series solutions, the other solutions may as well be written as series of generalized 
Laguerre polynomials, since we have a —zero or a negative integer in U{a, b, y). 

In addition to have changed the notations, we have also reordered the solutions of [1] 
so that the pairs are obtained by using rules r2 and ri in the same sequence as in section 
3. Thereupon we find that the integral relations, the conditions for having terminating 
series as well as the asymptotic behaviour of each solution are the same as those in the 
corresponding pairs of section 3. 

First pair : v — irj in {U^, C/°^) 

U^{z) = e*'^^^^+^''-(-^2/2) ^oo^^ b!f^\-2iuz)''U{n + 1 + 2ir], 2n + 2 + 2iri, -2iujz), 



U^{z) = e^'--^-^^-(S2/2) j^oc^^ ^(i) t/ (n + iry + f , 2n + 2 + 2iri, f ) . 



(54a) 



iuB^[n + 1] [n + 2 + i?7 - (^2/2)] 
2 [n + 1 + ^77] [n + i77 + (3/2)] ' 

/ . B2\ r]ijBi[{B2/2) - 1] 



/3« -53+ (n + l + ir]-^) (n + ir] + ^ 

(1) _ iuBi [n + 2iri] [n + {B2/2) + iri - 1] 
~ 2 [n + iri][n + iri- {1/2)] 



{n + irj) {n + 1 + ir)) 



(54b) 
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Integral relation (22): U[{B2/2) - i?7 - 1] > and ^{Bjz) > 0. (54c) 

r Eq. (51), if ir^^ 0,-1/2; 
Recurrence relations : < Eq. (52), if i?] = — 1/2; (54d) 

^ Eq. (53), ilir] = 0. 

Finite series, if (-B2/2) + irj = 1 - N =^ 0<n< N -1. (54e) 



The asymptotic behavior is given by equations (40a-b). 

Second pair : u = it] in {U^, U21,) or oc r2U^ and U2 oc r2U^- 
U^{z) = /(^)^i+ir/-(S2/2) ^00^^ h^^\-2iujzYU{n + 1 + 2ir], 2n + 2 + 2ir], -2iujz), 

Ul{z) = f{z)z-'^-^'''/^^ E~ €^ {-^y U(n + 2 + irj-f,2n + 2 + 2ir], ) , (55a) 



iuBijn + l\[n + iB2/2) + iri] 
2[n + l + i7]][n + iri + (3/2)] ' 

52\ / , . , 52\ , 77^^Si [(S2/2) - 1] 



(2) _ [n + 2^?7] [n - {B2/2) + ir] + l] 



2[n + i77][n + i77-(l/2)] 



Integral relation (27): 3?[(52/2) + ir/ - 1] < and ^{Bjz) < 0. (55c) 

r Eq. (51), if zr^T^ 0,-1/2; 
Recurrence relations : < Eq. (52), ifi77 = — 1/2; (55d) 

^ Eq. (53), if i77 = 0. 

Finite series, if (S2/2) -ir)^N + l 0<n<N -1. (55e) 
The asymptotic behavior is given by equations (42a-b). 



Third pair : u = S2/2 - 1 in (C/~, C/°^) or [71° oc riC/° and [7° oc riC/~ 

U^{z) = e*^^ E^=o &i^H-2icuz)'^?7 (n + ^ + ir/, 2n + ^2, -2icuz 
[/^(z) = e-^^i-^^ E^=o i'i'^ (f U{n + B2-I,2n + B2, f ) . 



(56a) 



a. 



(3) 



icc;5i[n+ 1] [n + (S2/2) - ir/] 



2[n+(52/2)][n+(S2/2) + (l/2)]' 
M3) r r^cui^i [(i?2/2) - 1] 

/3( ) = 53 + n(n + ^2 - 1) + ^ (^^/2) - 1] [n + (^2/2)] ' 

(3) _ icu^i [n + ^2 - 2] [n + {B2/2) - 1 + ir]] 
^" " 2 [n + (i?2/2) - 1] [n + (^2/2) - (3/2)] ' 



(56b) 



Integral relation (22): 3f?[(S2/2) - ^77 - 1] > and ^{Bi/z) > 0. (56c) 

r Eq. (51), if 52^ 1,2; 
Recurrence relations : < Eq. (52), if i?2 = 1; (56d) 

[ Eq. (53), if ^2 = 2. 

Finite series, if (^2/2) + irj = 1 - N ^ 0<n<N-l. (56e) 
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The asymptotic behavior is given by equations (44a-b). 



Fourth pair : u ^ 1 - B2/2 in {U^,U^^) or oc r2U^ and oc raC/l". 

U^{z) = f{z)z^-^^ En=o h^nK-'^i^zYU (fi + 2 - f + i7] , 2n + 4 - B2, -2iujz) , 
U2{z) = f{z)z-' E^=o b^r^^ (-f u{n + 3-B2,2n + 4-B2,-f), (57a) 

(4) _ icoB^ln + 1] [n + 2 - {B2/2) - irj] 

- 2[n + 2- {B2/2)] [n + (5/2) - {B2/2)] ' 

. - (n + + 2 - - (5Tb, 

(4) ^ w^i [n + 2 - ^2] + 1 - (^2/2) + ir/] 

- 2 [n + 1 - (^2/2)] [n + (1/2) - {B2/2)] ' 



Integral relation (27): 3?[(52/2) + ^r^ - 1] < and ^(^i/;^) < 0. (57c) 

r Eq. (51), if ^2 7^2,3; 
Recurrence relations : < Eq. (52), if B2 = 3; (57d) 

^ Eq. (53), if ^2 = 2. 

Finite series, if (S2/2) -ir)^N + l 0<n<N -1. (57e) 



The asymptotic behavior is given by equations (46a-b). 

We remark that, in addition to the three possible forms for the recurrence relations, 
the coefficients of the latter are fractional. Thence, these relations are not well defined 
when some denominator vanishes. Thus, if ir) — negative integer or half- integer < —1/2, 
the coefficients of the first and second pairs are not well defined, but we can form well- 
defined expressions by using the rule r^ which changes (77,0;) by (—77, —a;). Similarly, if 
for some value of B2 a denominator vanishes for the the third pair, we must consider the 
solutions of the fourth pair, and vice- versa. 

5. DCHE and GSWE: special cases and examples 

In this section we examine the two differential equations (lOa-b) which share the property 
of being particular cases of both the DCHE and the GSWE. In such equations, namely, 

+ [^0 + Oi cosh{Ku) + 62 cosh(2ft:'u)] Wi = 0, (58) 



2 



dv? 



+ 



^o + ^isinh(Kw) + e2COsh(2«;ii) W^2 = 0, (59) 



K is a given constant such that ku is real or pure imaginary and the 9i {6i) are con- 
stants. Thus there are only three parameters in each equation. If ku is pure imaginary, 
the first is the Whittaker-Hill equation (WHE) and if ku is real, the modified WHE [14] . 
In fact, Decarreau, Maroni and Robert [8] have already found that the WHE has that 
property, whereas we have found that these two equations are special cases of the GSWE 
[1]. Now we find some normal forms for the DCHE, from one of these we get the partic- 
ular equations written above and show that they also come from a GSWE with 1^2 = 1 
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and Bi — —Zq/2. Finally, we discuss solutions for problems obeying DCHE and GSWE, 
intending to decide on the best interpretation for each of the special equations (58-59). 

5.1. Normal forms for the DCHE 

Several normal forms for general Heun's equation and its confluent cases (except the 
triconfluent equation) were established by Lcmicux and Bosc [15]. Below, we give some 
forms suitable for DCHE when it is written as in equation (1). 
By performing the substitution 

U{z) = z-^^/2e^i/('^)F(z) (60a) 

in equation (1), we find for f{z) an algebraic normal form of the DCHE , namely. 



dz^ 



+ 



2 2?7CJ 1 / 
h — I S3 



B2 B2 



Bi _ B2 



Az* 



F = 0, (60b) 



z z^ \ 4 2 J z^ \ 2 
where, as before, Si 7^ 0, a; 7^ 0. The further transformations 

z = e^", F{z) = e^"/2py(M) ^ W{u) = z^^'-^^^''e-^'/'-^'^U(z), (61a) 
where A is a constant at our disposal, bring the equation to the hyperbolic normal form 
(PW 



dv? 
I{u) 

Si (1 



+ X'l{u)W = 0, 



(61b) 



s. 



— 2r)LJ 



2r)LV 
cosh(Aii) + 



sinh(Aw) + 



S2 



sinh(2Aii) + 



cosh(2Aii) + S3 - ^ (1 - S2)' 



to be used soon. The last form we shall need is obtained by the transformations 
z = p^, U{z) = p(i-2^^)/2e^i/(2''')G(p) ^ G{p) = z(2B2-i)/4g-Bi/(2.)^(^) 

in equation (1), and this affords another algebraic normal form given by 



(61c) 



(62a) 



d'^G 



dp^ 



+ 



.99 4 / B' B2 3 \ 



4Si _ S2 
2 



BI 



G' = 0.(62b) 



5.2. The common special cases 



Equation (61b-c) gives the relations among parameters which lead to equations (58-60) 
as special cases of the DCHE. For the WHEs we have 



-Si 



1 _ B2 
2 



^ + A2 [S3 - i (1 - B2f - 4770; cosh (Ali) + 2u;2cosh(2Aii)] = 0, 
and for the second equation 



= , 2r]u = Si 



1 _ 
^ 2 



du2 



+ A2 



S3 - i (1 - B2f - Ariujsm\i{\u) + 2^"^ cosh(2AM)l 1^2 = 



(63) 



(64) 



17 



Now we consider the particular GSWE 



z{z -^o)-^+ (^-y + + - '^Mz - zo) +uJ^z{z - zo)] [/ = 0, (65) 

\B\ — —zo/2, B2 — l\n equation (4)] 

which has only three constants, since Zq may be chosen at will, excepting zero. Then, the 
two special cases follow from the last equation by a change in the independent variable. 
For the WHEs we have 



z = Zq cosh^(crM/2), U{z) = Wi{u) ^ 

B3 + riLozo — ^uj'^zl — rjujzo cosh.{au) + ^uj'^Zq cosh.{2au) Wi — 0, 
and for the second equation 

z = f[i smh{au) + 1] , U{z) = W2{u) ^ 

B3 + rjUJZo — Ico'^Zq — irjLozo smh.{au) — kco^Zg cosh.{2au) W2 — 0. 



(66) 



(67) 



Note that the symbols Bi, rj and a; are denoting different objects depending on whether 
they appear in the DCHE or in the GSWE. 

In general, solutions for the WHE, when obtained from sohitions for the GSWE, are 
even or odd with respect to the variable u. For example: (i) wc have found pairs of even 
or odd solutions constituted by one solution in series of hyperbolic or trigonometric func- 
tions — Arscott's solutions [14, 16] — and another in series of Coulomb wave functions 
[1] , (ii) the Hylleraas and Jaffe type solutions to the GSWE [4] yield only even solutions 
for the WHE, but the transformation rule T2 given in equation (6b) of [1] generates new 
solutions whose limits are odd. 

5.3. The examples 

We write the one-dimensional time-independent Schrodinger equation for a particle with 
mass m and energy E as 

fp^ + [E-V{umu) = Q, u:=ax, ^ := (68) 

where a is a constant, x the spatial Cartesian coordinate, and ip must satisfy the regularity 
conditions 

lim ib = 0. (69) 

For the two following examples, we will see that it is convenient to interpret the WHE as 
a GSWE and the second special equation as a DCHE. 

First example: For the symmetric double-Morse potential of Zaslavskii and Ulyanov [17] 

B^ f 1 \ 

V^u) = — sinh^ u-B(s + -] cosh u, B > 0, (70) 

the Schrodinger equation is a modified WHE. If s is a non-negative integer or a half- 
integer, the potential (70) is a quasi-exactly solvable (QES) potential in the sense that 
one part of the energy spectrum stems from finite-series solutions [18]. If this WHE is 
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considered as a GSWE we find [1]: (i) even and odd quasi-polynomial solutions satisfying 
the regularity conditions, and (ii) even and odd infinite-series solutions which also satisfy 
the regularity conditions, provided that wc match solutions having the same series coeffi- 
cients and different radii of convergence. However, if the WHE is interpreted as a DCHE, 
only finite-series solutions satisfy the regularity conditions, as we will see next. 

For the sake of generality, we will get the solutions for the symmetric potential (70) 
as limits of the solutions for the asymmetric potential [17] 



V{u) = — (sinhii - — j - S (s + - j coshzi, > 0, C > 0) 



for which the Schrodinger equation becomes the DCHE 



+ 



CB 



1 



^2 



sinh u + B[s + -] cosh u — cosh(2M) + £ ^ + 



8 



8 



Comparing the above equation with equations (61a-c), we find 



A = 1 

and the parameters 
B 



z = e", ^{u) = W{u) = z 



B 



{B2-l)/2-Bi/{2z) 



U{z), 



(71) 



ij{u) = 0. (72) 



(73a) 



B 



C 1 



Si = -, ^2 = 1 + C - 2s, ^3 = £ + — + - sC, w = iry = -- - - - s. (73b) 

To form regular solutions we identify U [z] with the first and the third pairs of solutions 
given in section 3. The first pair yields 



^~(«) = g-f cosh«+(£+5)u ^(1) (^le'^) " , 

^^^{u) = e-f ^°^'^"+(^+^)"E~ o&^'^t^ (n - 2s,n + 1 - C - 2s, f e"") , 
with the following coefficients in the recurrence relations (29b) for h^^ 

a« = -(n + l), 

= -8- s{s + C) 



(74a) 



n[n 



C-2s) = P^'^ 



(74b) 



In 



— in -2s -I) 
4 



If s is a non- negative integer or half-integer (QES potential), we have 72s+i = and thus 
the series are finite with < n < 2s. Thence, the recurrence relations can be written as 



f Pq ao 

7i Pi «i 

72 ^2 

V 



\ 





Q!2 



723 P2s J 



( bo \ 



b2s-l 



/ bo \ 



£ 



b2s-l 
V b2s 



(75) 



This system of equations determines 2s -|- 1 different and real values for £ since, for a 
eigenvalue problem like this, that is, with a tridiagonal matrix, the following theorem 
holds (see [14], page 21): 'if aj, f3p 7^- are real and each product Oij'Jj+i is positive, then 
the roots corresponding to the equation (75) are all real and different'. In fact, this 
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remains valid even if both aj and 7^ are pure imaginaries, since we can take c„ = and 
put aj = —iaj, 13 = 7^ = 27^ in the recurrence relations for c„. The second condition 
of the theorem stands for the present case because 

q;,-7,+i = -{j + - 2s)SV4 > 0, 

where the inequality follows from the fact that < j < 2s — 1 for the elements of that 
matrix. Note that for these quasi-polynomial expansions we can select either the solution 
ip^{u) or ip^{^), since the hypergeometric functions in the latter reduce to a generalized 
Laguerre polynomial. The same is true for the solutions that result from the third pair 
of section 3, 

^^(u) = e-f (n-2s,n + l + C -2s, f e") , 

(76a) 

for which the coefficients of the recurrence relations for b!^^ are 

Pl^^ ^-E- s{s -C)-n{n + C- 2s) = '^f^ - (76b) 
if = ^{n-2s-\). 

The above pairs of solutions are related to one another by means of 

(i,'^{C-u),i,\{C-u)) ^ (f^{^-C--u),i,t{-C--u)) , (77) 

that is, these pairs remain invariant under the change {C,u) <-> {—C,—u), a property 
also present in Schrodinger equation (72). Note that we have obtained in [1] one pair of 
finite-series solutions for equation (72) by using the first pair of solutions given in section 
4.2; another pair might be obtained from that by using equation (77) or, alternatively, 
by using the third pair of solutions given in section 4.2. However, such solutions have 
recurrence relations with fractional coefficients and, for this reason, they are not well 
defined for certain integer values of the parameter C, as noted there. 

If we suppose that s is not a non-negative integer or half-integer, we can — following 
Leaver [4], section 8 — match the infinite-series solutions of each pair to get regular 
eigenfunctions which converge over the entire range of u, since the i/j'^ converge when 
exprt > and the the ■0° converge when expu < 00. In this case, the energy spectrum 
may be determined from the infinite continued fraction (30). Note also that the conditions 
to get from by an integral transformation are 3fJ(C) > and K(i?e~") > and, 
therefore, are assured. 

On the other hand, by using the second and fourth pairs of solutions given in section 3 
we may form two pairs of infinite-series solutions (even if the potential is QES) but these 
are not regular when u ^ —00 due to a multiplicative factor exp [—(B/2) sinhw]. 

For the symmetric case {C — 0) the two pairs degenerate to only one 

^O(k) = g-|coshn+s« ^00^^ (n-2s,n + l- 2s, f e-") (78a) 
= (5/2)2^6-1"°^^"-^'' E^=o ^i^^ (f " [see equation (47)] 



20 



and, in the recurrence relations (75), we have 



.(1) = -I 



a- = -(n + 1), ^-S-s'- n{n - 2s) = P^'^ - S, = ^(n - 2s - 1) (78b) 



52 



For finite-series (QES potential) the two solutions are convergent and regular for u G 
{—00, 00), but are neither even nor odd with respect to u. However, we can form even 
and odd eigenfunctions by taking the linear combinations 



2s 



n=0 



2s 



and choosing the constants Ai and A2 such that we have series of cosh [(n — s)u\ and 
sinh [(n — s)m], respectively. Note the absence of infinite-series solutions for this QES 
potential, contrary to the solutions resulting from the GSWE [1]. 

Second example. This example also deals with a QES potential, now giving a differential 
equation of the second type. In contrast with the the first example, we find that: (i) 
if the equation is treated as a DCHE, only infinite-series solutions satisfy the regularity 
conditions, (ii) if the equation is treated as a GSWE, there is no regular solutions. 
The potential is 



V{u) 



4 



/ 1\ 13 
sinh^u — (s + 7^ ) Bsmhu, B > 0, s = 0, -, 1, -, 



and thus the Schrodinger equation (68) reads 



+ 



E -I h \s^ - \B sinh u cosh(2M) 

8 V 2/ 8 



Comparing this equation with equations (61a) and (64), we get 
together with the following expressions for the parameters 



B B"^ 
B^ = --, B2 = 1 - 2s, B3 = S + — + s\ luj 
Z o 



B 



ir] 



s. 



(79) 



(80) 



(81a) 



(81b) 



Now, using the second and fourth pairs of sohitions given section 3, we obtain two pair 
of infinite-series solutions for '0(u). The first is 



^O(ti) = e-T Y^oo^^ h'^^^U (n+l,n+l- 2s, f e"") , 



(82a) 



with 



a«=n + l, p^^) ^g + s' + n{n-2s), 7!'^ = 



(1) 



-n. 



(82b) 



in the recurrence relations (29b). The second pair is given by 

(xi) = e-t cosh«+(i+s)« ^00^^ i^{4)u (n + l^n + 3 + 2s, f e") , 
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(83a) 



and has 



(83b) 



in the recurrence relations (29b). In these two pairs of nonterminating series we have to 
match the solutions in order to assure regularity over the entire interval for u. Now, the 
conditions to derive t/j^ from t/j^ by an integral transformation are 5R(2s + 1) > and 
3fJ(i?e~") > 0. Quasi-polynomial solutions result from the first and third pairs of section 
3 but they are not regular due to the presence of the factor exp [—{B/2) sinhu], as in the 
previous example. 

We have not found any regular solution by describing this problem by a GSWE. In 
effect, comparing equations (80) and (67), we find 

a — Zq — 1, z — - smhu+ -, cu — B, ir) — —s , — o -\ — t [s -\ — ] B. 



However the solutions for the GSWE have the factor exp {±iujz) = ±i?(— sinh-u + i)/2 
which diverges when m — > cxd or m — > — oo . As an illustration, we use the pair of solutions 
given in equation (42) of [1] and obtain 

~ oo 

^1 = g-(Bsinh«)/2 ^ ^(1) (5 sinh u - iByU (n - s, 2n + 1; S sinh u - iB) , 

n=0 

where the recurrence relations have the form given in equations (52), with 



iB 



2 



a- = -^(n + s + l), pi^)^-n'-S-^^ ^(i) ^ (^ _ s _ 1) . 



iB 

T 



Hence it follows that, if s is a non-negative integer, we have finite series and, if s is a 
half-integer, we have infinite series, but none satisfies the regularity conditions. 

We note that the potential (79) was obtained from the potential (70) by means of 
the change [s + (l/2)]Scoshw — >[s -|- (l/2)]Ssinhii, but it can also be derived from the 
potential 



\ TA • 1 2 Tr ■ , Vasmhu + V^ , , 

V \u) = V\ smh smh u H § 1 ^ \~^^ ^) 

cosh u 



(84) 



for which the Schrodinger equation is a GSWE [15]. A QES version of (84) is given in 
[19] and from that we can obtain (79) by taking V3 = V4 = 0. 

To finalize we mention two other problems which are reducible to a DCHE in one or 
another of the algebraic normal forms. They are related with the three-dimensional radial 
Schrodinger equation 



d^R 2 dR 
+ 



dr^ r dr 



+ 



2mE 

R — 0, t :— ^2 ' 



(85) 



for a particle with mass m and energy E. A normal form for this equation is 

d'^H{r) 



dr'^ 



+ 



H{r)^0, H(r) :^rR(r). 



(86) 
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Then, for the inverse fourth-power potential 



V{r) = Vir'^ + V2r-^ + V^r-^ + V^r 



-4 



(87) 



equation (86) assumes the normal form (60b) for the DCHE with z — r, whereas for the 
even-power potential 



it assumes the normal form (62b) with p — r. Solutions have been proposed for problems 
like these [20-23], but it would be intersting to study such problems from the viewpoint of 
the DCHE since the transformations rules and integral relations allow us to obtain new 
solutions from a known one. We could also check whether the solutions presented here 
are useful for that purpose. 

6. Concluding remarks 

We have found integral relations for solutions of the double-confluent Heun equation and 
combined them with transformation rules in order to obtain the group of solutions given in 
section 3. In section 4, the integral relations (22) and (27) have also been used to connect 
expansions in series of Coulomb wave functions. The solutions have been displayed in pairs 
{U°°{z),U^{z)), where U°°{z) and Uf{z) converge for \z\ > and \z\ < oo, respectively. 
The solutions in each pair have the same series coefficients just because U^{z) comes from 
U°°{z) by an integral relation, whose validity conditions (18a-b) were satisfied due to the 
choice of an integration contour that excluded the point t = 0, where the solution [/?° {t) 
does not converge. 

Comparing each pair of solutions given in section 3 with the corresponding pair given 
in section 4.2, we find (i) the same integral relations between the solutions, (ii) the same 
asymptotic behaviour for the solutions, (iii) the same conditions for quasi-polynomial 
solutions, and (iv) the possibility of writing one solution as a generalized Laguerre poly- 
nomial. The last property may be important if we need to normalize solutions. Note, 
however, that in section 3 there is only one form for the recurrence relations, whereas in 
section 4.2 there are three possible forms with fractional coefficients which are not well 
defined when a denominator vanishes. 

In section 5 we have given normal forms for the DCHE and analysed equations (10a- 
b) that have the common property of being particular cases of both the DCHE and the 
GSWE. We have as well looked for solutions to the Schrodinger equation for two quasi- 
exactly solvable hyperbolic potentials. For the symmetric Zaslavskii-Ulyanov potential we 
have a modified WHE, for which we have established regular quasi-polynomial and infinite- 
series solutions by considering this WHE as a GSWE [1], but only quasi-polynomial 
solutions by treating it as a DCHE. For the potential (79), in the second example, the 
Schrodinger equation is an equation of the second special type, for which we have found 
regular infinite-series solutions by interpreting the equation as a DCHE; however we have 
not found any regular solution by considering it as a GSWE. For infinite-series solutions, 
which have been formed by joining solutions belonging to a same pair, the energy spectra 
may be computed from infinite continued fractions, for instance. 

The results of the preceding paragraph suggest that we must interpret the WHE 
(10a) as a GSWE, and equation (10b) as a DCHE. However, it is necessary to consider 
other problems such as the time dependence of the Dirac equation in radiation-dominated 
Friedmann-Robertson- Walker spacetimes (see section 2.2.1 of [1]). The differential equa- 
tions for these problems have no free parameters, which implies that we have to deal 



V{r) = Fir^ + V2r-^ + V^r'"^ + V^r 



-6 



(88) 
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with double-sided series possessing a phase parameter, as those of section 4.1 for equation 
(10b) or the solutions of section 2.2 of [1] for the WHE. 

In section 5 we have also discussed solutions to the Schrodinger equation with an 
asymmetric Zaslavskii-Ulyanov potential, but now using the new solutions found in section 
3. The regular solutions are quasi-polynomial and do not exclude integral values for the 
parameter C, contrary to the solutions constructed on the basis of the expansions given in 
section 4.2 [1]. Despite this advantage, the obtainment of regular infinite-series solutions 
for this problem remains unsolved. 

Finally, we have called attention to some singular radial potentials for which the 
Schrodinger equation leads to DCHEs, but we have not tried to solve these equations. 
Another related problem concerns the solutions of the Schrodinger equation for the QES 
potentials derived recently by Bagchi and Ganguly [24], since for their hyperbolic poten- 
tials we will find modified WHEs, and for trigonometric potentials, WHEs. Considering 
these equations as particular cases of the GSWE, we may expect to find quasi-polynomial 
and infinite-series solutions as in the first example of section 5. 

In the Appendix we have rewritten some integrals in a form suitable for use in section 
3 and 4. Thus, excepting a correction to a misprint in a table of integrals, this appendix 
does not contain anything original. 



Appendix. Integrals used in sections 3 and 4 

The first equation is an integral representation for the irregular confiuent hypergeometric 
function U{a,b,z) given by [12] 

e-y\t - ly-H^-^-^dt = r{a)e-yU{a, y), [Ka > 0, % > 0]. (Al) 

The two following integrals are usually given in terms of irregular Whittaker functions 
W'^ft,/j(^)- By convenience, we have reexpressed them in terms of irregular confluent hy- 
pergeometric functions by using the relation [13] 

= e-'/Y^^'/'^U (^i - + 2/. + 1, y) . 
The first of these integrals is 

e-v{y - i)M-i[/ Q _ ^ _ A, 1 - 2A, ay) dy 

= r{ii)e-''a-^U Q - K - A, 1 - 2A - /i, a) , [% > 0, 3?a > 0], (A2) 

which results from [25] 

e-y/^(y - l)^^~'y^-.WU<^y)dy = r(/.)e-»/^a-^/^iy,_|, ,+H(a), [3?// > 0, 5Ra > 0]. 

Note a misprint on page 867 of [26] where we have VrK-(^/2), x-(n/2){o,) on the right-hand 
side of the above integral. The other integral is 

e-«^(y - ir-Y^^-'-'^U Q + A - 2A + 1, ay) dy 

= r(/i)e-"?7 Q + /i - K + A, 1 + 2A, a) , [3fJ/x > 0, 3?a > 0], (A3) 
which is equivalent to [26] 

e-«^/^(y - ir-Y-'''~'W^Aay)dy = r(//)e-«/2W^«_^, ,(a), [^pt > 0, 5fta > 0]. 
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